0. Código

>	restart; with(plots):

Warning, the name changecoords has been redefined

>	scan   := 'scaling = constrained, axes = none':

>	r10    := proc()  2*(evalf(rand()/10^11)-5);      end proc: rr     := proc()  2+ 4*(evalf(rand()/10^12));     end proc: shade  := proc()  local s;                   s := evalf(rand()/10^12,2)/4 + .3;                   COLOR(RGB, s,s,.9);           end proc:

>	CP := (h,k,r) -> plottools[circle]([h,k], r, thickness = 2, color = shade()):

  1. O Número

Como quadrados, mas triângulos diferentes, todos os círculos são similares aos outros. Cada círculo é apenas uma versão menor ou maiorde todos os outros círculos.

>	display( seq(   CP( r10(), r10(), rr() ), k = 1..10), scan);

>	display( seq(   CP( r10(), r10(), rr() ), k = 1..10), scan);

>	display( seq(   CP( r10(), r10(), rr() ), k = 1..10), scan);

O raio da circunferência para o diâmetro é o mesmo para todos esses círculos. Aquele raio é . Suponha que nós medíssemos vários objetos - ambos ao redor e de lado a lado, e pegássemos seus raios.... esse é um bom projeto ... latas de lixo, latas de soda, lâmpadas, anfiteatros....qualquer coisa redonda poderia ser medida. As unidades não importariam se elas fossem as mesmas para ambas as medidas.

Aqui estão algumas amostras.

>	C := 8; d := 5/2; C/d : % = evalf(%, 5);

>	C := 17; d := 16/3; C/d : % = evalf(%, 5);

>	C := 69;  d := 22; C/d : % = evalf(%, 5);

>	C := 300; d := 95; C/d : % = evalf(%, 5);

>	C := 472; d := 150; C/d : % = evalf(%, 5);

  2. Expansão Decimal de

Existem lendas que dizem que um matemático dedicou toda a sua vida calculando os dígitos de . É um conto interessante. Nos anos 60, durante a era das regras, um estudante, professor, ou engenheiro poderia usar frequentemente 3.14 como uma aproximação para . Embora computadores existissem e muitos dígitos já fossem conhecidos nessa época, eles não eram acessíveis a uma grande quantidade de pessoas. No meio dos anos 70, as calculadoras se tornaram disponíveis por vários dólares que poderiam calcular    em dez posições decimais. Mesmo hoje, a maioria das calculadoras sãolimitadas a 15 ou 16 posições decimais.

>	Pi = evalf(Pi, 16);

No final do século 16,  foi conhecido em 30 posições decimais. No final do século 18, ele foi conhecido em 140 posições. E no final do século 19, ele foi conhecido em 500 posições precisas.

>	Pi = evalf(Pi, 30);

>	Pi = evalf(Pi, 140);

>	Pi = evalf(Pi, 500);

Mesmo com o advento dos computadores pessoais no início dos anos 80 era tão automático que poderíamos calcular  em várias posições decimais - porque o computador armazena os dados em bytes. Uma fração 32 byte tem essa resolucção de mais ou menos 10 posições decimais.

>	1/2^32; evalf(%);

Hoje, com o Maple, você pode ver o que não se pensava a algumas décadas atrás. O Maple inclui algumas programações que transcendem os limites de armazenameno em bytes para calcular um número quase ilimitado de posições decimais. Simplesmente escolha quantas posições decimais de ¹ você qier, e como se ele fosse um gênio, isso se tornará realidade . 100 posições decimais, 1,000 posições, 2,500 posições.

Seção Opcional (Utiliza Pré-Cálculo)

Se nós considerarmos as posições decimais conhecidas de    como uma medida de tecnologia avançada, parece ter crescido exponencialmente. Vamos ver dois pontos : (ano 1600, 30 posições) e (ano 1800, 140 posições) para extrapolar uma função.

>	C := 'C': t := 't':  k := 'k': f := t-> C*exp(k*t); EQ1 :=      f(1600) = 30; EQ2 :=      f(1800) = 140; solve(EQ1, C) = solve(EQ2, C);

>	k := solve(%,k); EQ1; C := solve(%,C);

Aqui está a função que deve fornecer o número de posições decimais de  baseada no ano.  

>	f(t); simplify(f(t),exp);

De acordo com essa fórmula, vamos ver quantas posições decimais deveriam existir nos anos 2000 e ano 0.

>	floor(evalf( f(2000) ));

>	floor(evalf( f(0) ));

Parece que os antigos estavam à frente da tabela! Ou os anos escuros nos colocam atrás na tabela.

>	floor(evalf( f(1200) ));

  3. A História de

O número , ou alguma variação dele, existe por no mínimo vários milênios. Quem sabe se o inventor da roda pensava sobre o raio da distância que a roda viajaria  comparad com a largura da roda - a qual é . No entanto, o valor verdadeiro de  não chega a ser conhecido imediatemante.

      = 3?

O número  tem um significado especial na matemática e na história. Vários ancestrais tais como os Babilônios usaram o valor de 3 para .
How accurate is this approximation?

>	3/Pi; % = evalf(%,10);

Está 95% preciso. Se um velho arquiteto de um coliseu circular acreditou que o raio fosse 300 pés, de quantos pés de cerca ele iria precisar?

>	2*300*3; 2*300*Pi; difference = evalf( % - %%);

Simplesmente por caus ada imprecisão desse ou daquele " ", o arquiteto poderia precisar de 85 pés de cerca adicionais. Ele deveria ter ficado maravilhado.

      = 3 1/7?

Um "novo e melhorado" , foi encontrado pelos ancestrais Egípcios. Eles usaram o valor de 3 1/7, o qual é 22/7.
  

>	EPi := 22/7;

>	EPi/Pi;   % = evalf(%,10);

Essa aproximação é melhor. É apenas  4/100 de um percento! Se o arquiteto usou esse valor em seu coliseu, o quanto ele estaria distante?

>	2*300*EPi; 2*300*Pi; difference = evalf( % - %%);

Mais ou menos  3/4 de um pé. Isso não é tão ruim. E se ele quisesse calcular a distância que a terra viaja, usando o fato de que o sol viaja em uma órbita circular sobre o sol (atualmente eliptica) e o raio é 10,000 milhas.

>	2*10000*EPi; 2*10000*Pi; difference = evalf( % - %%);

 

     = 3 17/120?

Um Romano usou um valor melhor, 3 17/120, o qual é 377/120.
  

>	RPi := 3 + 17/120;

>	RPi/Pi;   % = evalf(%,10);

Essa aproximação é melhor. É apenas 1/425 de 1 percento! E se o arquiteto usasse esse valor no seu coliseu, o quanto ele estaria distante?

>	2*300*RPi; 2*300*Pi; difference = evalf( % - %%);

Mais ou menos 1/20 de um pé. E se ele quisesse calcular a distância que a terra viaja, usando o fato de que o sol viaja em uma órbita circular sobre o sol (atualmente eliptica) e o raio é 10,000 milhas.

>	2*10000*RPi; 2*10000*Pi; difference = evalf( % - %%);

 

     = sqrt(98694)/100?

Os Hindus usaram essa aproximação.
  

>	HPi := sqrt(98694)/100;

>	HPi/Pi;   % = evalf(%,10); 100*(1-evalf(HPi/Pi,10));

Essa aproximação difere mais ou menos 1/1000 de 1 percento! E se o arquiteto usasse esse valor no seu coliseu, o quanto ele estaria distante?

>	2*300*HPi; 2*300*Pi; difference = evalf( % - %%);

Mais ou menos 1/5 de um pé. E se ele quisesse calcular a distância que a terra viaja, usando o fato de que o sol viaja em uma órbita circular sobre o sol (atualmente eliptica) e o raio é 10,000 milhas.

>	2*10000*HPi; 2*10000*Pi; difference = evalf( % - %%);

  4. Outra Aproximação Fracional de

Existem mais aproximações fracionais de   . Vamos ver uma série de frações precisas.

>	FracApprox := 22/7; x := evalf( %); `Percentage Error` = 100* abs( evalf( (Pi-evalf(FracApprox))/Pi, 50) );

>	FracApprox := 333/106; x := evalf( %); `Percentage Error` = 100* abs( evalf( (Pi-evalf(FracApprox))/Pi, 50) );

>	FracApprox := 355/113; x := evalf( %); `Percentage Error` = 100* abs( evalf( (Pi-evalf(FracApprox))/Pi, 50) );

>	FracApprox := 103993/33102; x := evalf( %); `Percentage Error` = 100* abs( evalf( (Pi-evalf(FracApprox))/Pi, 50) );

>	FracApprox := 104348/33215; x := evalf( %); `Percentage Error` = 100* abs( evalf( (Pi-evalf(FracApprox))/Pi, 50) );

>	FracApprox := 208341/66317; x := evalf( %); `Percentage Error` = 100* abs( evalf( (Pi-evalf(FracApprox))/Pi, 50) );

>	FracApprox := 312689/99532; x := evalf( %); `Percentage Error` = 100* abs( evalf( (Pi-evalf(FracApprox))/Pi, 50) );

>	FracApprox := 833719/265381; x := evalf( %); `Percentage Error` = 100* abs( evalf( (Pi-evalf(FracApprox))/Pi, 50) );

>	FracApprox := 1146408/364913; x := evalf( %); `Percentage Error` = 100* abs( evalf( (Pi-evalf(FracApprox))/Pi, 50) );

>	FracApprox := 4272943/1360120; x := evalf( %); `Percentage Error` = 100* abs( evalf( (Pi-evalf(FracApprox))/Pi, 50) );

>

>

Assim, usando frações "modestas", nós podemos obter algumas muito precisas!